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Abstract 

We solve the problem of determining airfoils that approximate, in a least square 
sense, given surface pressure distributions in transonic flight regimes. The flow is 
modeled by means of the Euler equations and the solution procedure is an adjoint- 
based minimization algorithm that makes use of the inverse Theodorsen transform in 
order to parameterize the airfoil. Fast convergence to the optimal solution is obtained 
by means of the pseudo-time method. Results are obtained using three different pres- 
sure distributions for several free stream conditions. The airfoils obtained have given 
a trailing edge angle. 
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1 Introduction 


State of the art computational techniques allow the analysis of complex two-dimensional 
or three-dimensional configurations of aerodynamic interest, using models of increasing 
complexity. In the design of aerodynamic components however, simpler physical models 
have been used, in order to obtain analytical solutions or tractable computational problems. 
In some cases, the use of simpler models does not impair the validity of the results obtained. 
For example, transonic shock-free airfoils were designed using the full potential equation 
written on the hodograph plane: since shock-free solutions were sought, it was not necessary 
to use the Euler equations. 

This simplification is generally not valid when shocks, even weak shocks, are present in 
the flow field. In fact, although entropy production by weak shocks could be small, at shock 
points of high curvature, entropy gradients could be high as pointed out by Nieuwland & 
Spee [18]. Another reason to avoid the potential approximation in the presence of shocks 
is the possibility of multiple solutions as shown by Steinhoff Sz Jameson [24]. Therefore, 
since shock-free solutions are at least unlikely when designing airfoils in the presence of 
constraints of various kinds, the use of the Euler equations becomes necessary. 

Much effort in transonic aerodynamics has been devoted to the design of shock-free air- 
foils, even though according to an early theoretical result due to Morawetz [16], shock-free 
solutions of the full potential equations are isolated points. However, subsequent experi- 
mental investigation conducted by Pearcy [20], and by Nieuwland & Spee [17] showed that 
shock-free flows around airfoils with no or weak shocks can be experimentally realized and 
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are stable with respect to unsteady disturbances or variations of M <*,. The design methods 
employed at that time were based on the transformation of a known solution in the hodo- 
graph plane to the physical plane, and iterating on the input hodograph solution based on 
to the result obtained in the physical plane. For a review of these methods see the article by 
Boerstoel [4]. A common drawback of these approaches is that they have very little control 
on the resulting pressure distribution over the airfoil. Volpe [27] tried to overcome this 
problem by taking pressure distribution as an input design requirement. However, for the 
well posedness of the problem, special restrictions had to be applied on the possible choices. 
In addition, the flows obtained with such airfoils showed secondary shocks embedded in the 
flow field, that caused losses even at design conditions and boundary layer separation at 
off-design conditions, as pointed out by Labrujere h Slooff [15]. 

More recently, advances in design methodology have come about through shape optimiz- 
ation. Shape optimization is based on the minimization of an objective function, defined 
on the boundary, with respect to the variation of the boundary itself. This approach allows 
for example multi-point design and the easy introduction of design constraints, in addition 
to control of the pressure distribution on the body. The first attempts of Pironneau [21] [22] 
and of Glowinski h Pironneau[7] to solve such problems, were based on control theory. 
With this approach the minimum of the functional is determined with a descent algorithm 
based on the gradient of the functional with respect to the controls. The gradient of the 
functional is determined solving the adjoint of the equations governing the flow. No much 
progress was made at the time because no code for the computation of the flow field and of 
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the adjoint equations was available. Recently this method was reintroduced by Jameson [13] 
and computationally demonstrated its inherent vantages. In fact, with other minimization 
strategies for wing design, such as that first proposed by Hicks & Henne [9], the gradient of 
the functional is computed by finite differences, perturbing the controls one by one. This 
requires a complete flow-field solution for each gradient component to be computed. Even 
though optimal designs have been obtained with this method for the potential equations, 
see Reuther, van Dam k Hicks [23], the method becomes too computationally expensive 
for the Euler and Navier-Stokes equations. In contrast, with the method first proposed by 
Pironneau [21], it is necessary to compute the flow field only once for each gradient eval- 
uation. This has led to several applications for both internal and external flows governed 
by Euler equations - see Beux k Dervieux [2], Jameson k Reuther [14], Iollo k Salas [11] 
- for Navier-Stokes equations, see Cabuk k Modi [5]. In the latter work, however, the cost 
of computations was so high that they could only afford to compute about 10 flow-field 
solutions on a relatively coarse grid. 

A further decrease in the computational effort is possible by applying the theory proposed 
by Ta’asan [25], namely the pseudo-time method, which is based on the following observa- 
tion. Gradient-based methods (including adjoint formulations) can be viewed as marching 
along the intersection of the hypersurfaces representing the solutions of the flow-field and 
of the adjoint equations. This is an expensive process since each step requires the solution 
of two systems of PDEs. The idea of the pseudo-time method is to perform the march- 
ing while satisfying the optimality conditions on the boundary, but without satisfying the 
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flow-field and adjoint equations in the interior domain. The cost of such an iteration per 
step is significantly smaller than that of gradient-based methods. Its convergence has been 
shown by Ta’asan to be independent of the number of design variables. This method was 
applied by Iollo, Kuruvila & Ta’asan [10] to a few model problems governed by the Euler 
equations, which require the solution of the optimization problem in a small vicinity of the 
boundary. The results were compatible with the theoretical predictions. 

In the present paper we use the above method for transonic airfoil design. In particular we 
are concerned with the problem of finding families of wing sections that have approximately 
the same pressure distribution at different flight conditions. An example of such airfoils 
was given by Chin [6]. In this paper we show a few families corresponding to different 
pressure distributions, with prescribed trailing edge angle. Such families are designed by 
means of the inverse Theodorsen transform, that allows the optimization with an increasing 
number of harmonics defining the airfoil geometry. Compared to previous work, we are able 
to design for given pressure distribution using the Euler equations, even in the presence of 
flows with sharply captured shocks. The total computational cost to reach an optimum is of 
the order of the cost of a few flow-field analyses, allowing a cheap mapping of performance 
with several target pressure distributions. 

2 Statement of the Problem 

Consider a plane transonic flow around an airfoil as in fig. 1. Let the upstream Mach number 
be Moo and take the unperturbed stream pressure and density as reference. Disregarding 
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Figure 1: Physical plane. The portion of the plane included in H, and external to F is 
denoted by ft. 

viscous effects, the flow is governed by the Euler equations written in the following form 

Uf + F x + G y = 0 ( 1 ) 

where 

p pu pv 

pu p-\- pu 2 puv 

U= F= G = 

pv puv p + pv 2 

pe u(pe-Fp) v(pe + p) 

Here u, v, p, e and p are the velocity components, density, specific total energy and pressure 

respectively. Assuming a perfect gas law, we also have p = Kp( 2e — u 2 — v 2 ) with k = 

(7 — l)/2 and 7 the ratio of specific heats. The Euler equations define the state of the 
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Figure 2: a )z plane; b)to plane; c)w" plane. 

system, and therefore are named state equations. 

We assume M 0 0 < 1 and that the external boundary S is far enough so that the flow is 
unperturbed there. If the flow is incoming with respect to the normal to E, total pressure, 
entropy and the ratio v/u are assigned, whereas if the flow is outgoing static pressure is 
prescribed. The boundary condition on the airfoil is p V • n = 0, where n is the normal to 
F oriented as in fig. 1 and V = («,»)• 

Consider a real valued functional £(F,X(F)) where X(F) is a solution of the Euler equa- 
tions with boundary conditions. The optimization problem consists of minimizing the func- 
tional 8 over all the admissible shapes of the boundary F. 

The functional 8 is chosen according to some design criteria. Although the formulation 
we propose is general, here we will limit ourself to the case in which 

£(T) = lf r \p(T)-f?ds ( 2 ) 

where p* is a given pressure distribution and s is the curvilinear coordinate on F. The 
solution of such problem will allow the design of minimum-wave-drag profiles in transonic 
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Figure 3: Computational mesh. 


flows. Constraints on aerodynamic characteristics of the airfoil such as pitching moment or 
lift coefficient can be added to E in the form of quadratic penalty functions. 


3 Analysis 


In this section we will derive the necessary conditions for the minimum of £(T) subject to 
the Euler equations. The derivation is similar to that presented for internal flows by Iollo 
& Salas [11]. Consider the augmented functional 


£(U, r, A, p) = € + f *A(AU a + BU y )dQ + / ppV 

Jtl jT 


• nds 


(3) 


where A = dF/dU and B = dG/dU. The vector A (x,y) = *(Ai, A2, A 3 , A4) and the scalar 
are the continuous equivalents of the Lagrange multipliers. We will show that the 
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unconstrained minimum of eq. 3 with respect to [/, T, A and p is the same as that of £(T). 

Suppose each point of the airfoil can be displaced a length £n(s) in the direction of the 
normal n to the profile in a way that will be described in the next section. Corresponding 
to this change, U(x,y) is increased by a function sU(x,y); A (x,y) by eA(x,y) and p(s) by 
ep(s). Calculating the variation of the functional C with respect to the variation of U , A, 
y and f respectively, and disregarding higher order terms we obtain 


where 


and 


SCu 


jf J|(p( r) - f)Uds + J 'A {An x + Bn y )Uds + 
Jj?A x A + 'AyBpdCl + J /in • ^ Uds 





and 


dpV 

dU 


( \ 

0 10 0 

v 0 0 1 0 } 


(4) 


SC A = 

[ *A (AU X + BU y )da 

(5) 

II 

a. 

‘-o 

j) ppY • nds 

(6) 

SCr = 

k - p") s ds + L <A ( au * + bu v ) « ds + 



J p ^ ^ « n nds J p pY • n ds 

(7) 


where n is the variation of the normal to the airfoil boundary T. 

At the minimum of the functional, for all the possible choices of the functions [/, A, p and 
h, the following must be true 


SjCjj = SjC\ = — SjC p = 0 


(8) 
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Therefore, we have 


SCa = 0 <£> AU X -f BU y = 0onfi 

and 

SC M = 0 ^ pV • n = 0 on T 

which are the Euler equations and boundary conditions. Furthermore we obtain the follow- 
ing adjoint or costate equations 



( 9 ) 


( 10 ) 


where 

p = ~ [^1 + u\ 2 + v\ 3 + (je - kV 2 )\ 4 (11) 

From eq. 10 the boundary condition for A on F are found. Given U, the set of costate 
eqs.(9-ll) determine A in 0 and p on F uniquely. 

It is seen that if we could solve for h from the system of eqs. 8, then we would also 
determine the minimum of S. Unfortunately, eqs. 8 cannot be easily solved for h. Nev- 
ertheless, using an optimization algorithm it will be possible to determine the shape of T 
that minimizes C and S. In fact, knowing U and A for a given shape of the airfoil, we can 
calculate from eq.(7) how to choose h so that SC < 0, and iterate until the minimum is 
reached. Two algorithms of this kind will be discussed in Section 5. 
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In case of shocks occuring in the flow field, we split the domain of integration by means 
of a curve T that coincides with the shocks where they exist. Then we follow the same 
derivation presented so far on each of the snb-domains, regarding T as a boundary; the 
situation is similar to that which occurs for internal two-dimensional flows. The resulting 
extra condition for A on the shock is A = 0. It should be noted that if the shocks are not 
treated properly, the problem of solving the costate equations with boundary conditions is 
not well-posed, as shown by Iollo, Salas &; Ta’asan [12]. 


4 Airfoil Design by Inverse Theodorsen Transform 


By means of a conformal mapping, it is possible to transform a large class of shapes of 
aerodynamic interest to the unit circle. In particular, using the Theodorsen transform it is 
possible to transform any given profile to the unit circle. As an example, we will transform 
the NACA 0012 profile of fig. 2a. First the profile is represented in the complex plane z, 
with the trailing edge corresponding to the point (1,0). Then the plane z is transformed 
into the plane w (fig. 2b) using the relation 

z-1 

W = r 

Z + 1 

The angle in the origin is eliminated by the transformation 



where 


/ T 

w = w 


IT 


T = 


2tt 


a 


( 13 ) 
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where a is the angle at the trailing edge. Finally, the plane w' is transformed into the 
quasi-circle of fig. 2c by 


w 


n 


l-tc' 
w' + 1 


(14) 


If we denote by f the unit circle plane, the transform is completed by 


w " = 


(15) 


where the coefficients a n and b n are usually determined by means of the Fast Fourier Trans- 
form. 

In the design of an airfoil, we proceed in the opposite way. Starting from the unit circle 
plane, the coefficients a n and b n are given so that an airfoil is obtained on the £ plane. Let 
us define £ = r e zv> and w" = pe l 9 . In order to design an airfoil with a given trailing edge 
angle, it is necessary that the figure resulting on the plane w" passes through the point 
(1,0). In addition it is convenient that for = 0 we have 0 = 0. It is easily seen that for 
(p = 0, the radius r is not affected by the value of 6 0? while 0 = Elo &n* Hence, taking 
bo = — Eli we bave <p = 0 for 0 = 0. Finally, since the value of 0 is not affected by 
the coefficient ao, we can take do = — Eli a n ' m order to have p = 1 for ip = 0. 

The coefficients must be evaluated accurately to obtain reasonable shapes in the 2 plane. 
For example, in the case of non lifting profiles b n — 0 for every n. On the other hand, the 
summation Elo n ^ controls the camber, as can be seen studying the transform in the 
vicinity of (p — 0. 

Given a set of coefficients a n and b n we are now interested in calculating the variation of 
the tangent and of the normal when we increment one of the coefficients by a small quantity, 
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say Scii. Given two points 2: and z' on the airfoil, the tangent to the profile is defined by 

z' — z , 

lim r— : ( 16 ) 

| z* - Z\ V 7 

Let us consider the corresponding points on the unit circle, respectively f' and f , we have 


lim ^ |r = e ‘ (*H-»/2) 

«'-* If - £1 


Therefore we obtain 


lim ^ e i(v+ ' /2) ( 18 ) 

1 2?' - 2f| df V 7 

where H is the transformation from f to 2?. For example, the variation of the normal with 


respect to a,- is given by 


~ 5 ( dH \ . 

n = e if ) Sai 


Often airfoils are subject to geometric constraints such as fixed volume, maximum or min- 
imum thickness, or given trailing edge angle. These constraints can be easily included in 
the mapping formulation. For example, if r is kept constant, then the trailing edge angle 
is constant. 

On the other hand, the volume included by an airfoil is 


-*(/ 


and hence it is possible to compute the variation of the volume with respect to the coefficients 
of the transformation. Therefore, when the volume is kept fixed, the equation 


eS^+eS^q 




n=0 uu n 
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represents the additional constraint in the minimization. In particular, we can take for 


example 


E£U,„* (dV/da n ) 5a + £"=0 ( 9V/8K ) Sb 


N 


8 a; = 


( 22 ) 


dV/dat 

and substitute Sai in the expression for the descent direction, in order to obtain a constrained 
minimization. 

Taking the coefficients a n and b n to be the controls, we were able to perform a minim- 
ization with an increasing number of controls involved. In fact, in the first optimization 
iteration, the rough shape of the airfoil must be defined, therefore it is needless to work 
with many Fourier coefficients at this point. As the rate of convergence to the minimum de- 
creases, higher harmonics are added until the convergence rate increases and the functional 
is minimized. 


5 Numerical Solution 

In this section we will discuss the numerical solution of the flow-field equations and of the 
adjoint equations. In addition, two solution strategies for the minimization problem of 
Section 2 will be presented. 

The Euler equations are solved by means of an explicit, finite- volume scheme based on the 
flux-difference formulation proposed by Pandolfi [19], and the steady solution is obtained as 
an asymptotic limit of a time evolution. The conservative variables U are computed at the 
cell centers, and the fluxes F and G are evaluated at the cell interfaces with an approximate 
Riemann solver. Second-order accuracy is achieved using a method in the spirit of the 
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1 0.0500 

2 0.1000 

3 0.1500 

4 0.2000 

5 0.2500 

6 0.3000 

7 0.3500 

8 0.4000 

9 0.4500 

10 0.5000 

11 0.5500 

12 0.6000 

13 0.6500 

14 0.7000 

15 0.7500 

16 0.8000 

17 0.8500 

18 0.9000 

19 0.9500 

20 1.0000 

Figure 4: Iso-Mach lines for a NACA 0012 at 1 degree of incidence and M 0 0 = 0.85. 
C L = 0.334, C D = 0.0568. 

essentially non-oscillatory scheme of Harten, Engquist & Chakravarthy [8]. With such 
an approach the flow-field values at the cell interfaces, used as initial conditions for the 
Riemann problem, are reconstructed by means of a linear interpolation. The occurrence of 
spurious oscillations is prevented using a minmod limiter. The amplitude of the integration 
step is chosen in accordance with the Courant-Friedrichs-Lewy (CFL) condition. 

Compared to other solvers, this one allows a sharper capturing of shocks and contact 
discontinuities, so that wave-drag computations are more reliable. It should be emphasized 
that in other works, such that of Jameson h Reuther [14], no special treatment of the 
costate equations is made at the shock. This is allowed in their formulation because shocks 
are smeared over several grid points due to artificial viscosity. However, because of the 
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Figure 5: Pressure versus chord fraction for Korn airfoil and optimal airfoil at = 0.75. 
The unperturbed flow static pressure is taken as reference. 

smearing of the shock, the computed gradient might not be reliable close to the minimum, 
causing a drastic reduction of the convergence rate. 

The costate equations have no conservative form, and the numerical solution is obtained 
with a finite- difference scheme. We introduce a set of curvilinear coordinates (p(x, y) and 
'ip(x^y). The costate equations are then rewritten as follows 

*AA V + l BK^ = 0 (23) 

where A — Aip x + B<p y and B = Atf} x + Bijj y . The transformations c p and 0 are defined 
as (xi,y m ) l and (xi, y m ) A- m, where l and m are the indices that characterize each 
grid point. The above equations are linear and as such are their boundary conditions, see 
Iollo h Salas [11]. Suppose that locally we separate the variables through the following 
approximation: 

A (y>, 0, t) = A \(p, t) + A"(0, t) (24) 
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Figure 6: Pressure versus chord fraction for Korn airfoil and optimal airfoil at M^ = 0.8. 
The unperturbed flow static pressure is taken as reference. 

We are left with two one-dimensional subproblems that are solved as the asymptotic limit 
of a time-dependent technique, as explained by Iollo & Salas [11]. 

The computational grid is the same for Euler and adjoint equations and is obtained by 
means of the transformation that we use to generate the profile. All the computations have 
been performed on an O-grid with 120 tangential and 30 radial points, as partially shown 
in fig. 3. The outer boundary is at 10 chords from the airfoil. 

In Section 3 we have introduced the basic idea for the minimization of the functional £. 
From eq. 7 and eq. 19, once the solution of Euler and adjoint equations is computed for the 
present geometry, it is possible to calculate the gradient of £ with respect to the coefficients 
a n and b n . In particular, disregarding higher order terms, we have 

N 

sc = y ; 

71=1 


dC - dc ' 
rs OCLn T rj, bb n 

O&n 
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Figure 7: Korn airfoil and optimal designs at Moo = 0.75 and M^ = 0.8. 


and taking 

Sa n = —vdC/da n (26) 

and 

8b n — —vdC/db n (27) 

it is 8C < 0 if v is small enough. A first minimization algorithm is therefore: 

1. start with a set of a n and b n \ 

2. enforce 8C& = 0 and <5X M = 0 by finding a U that satisfies the steady state Euler 
equations and boundary conditions; 

3. enforce SCjj = 0 by finding a A that satisfies the costate equations and boundary 
conditions; 
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Figure 8: Pressure distribution over a Korn airfoil at M 0 0 — 0.8. 

4. calculate dCjda n and dC/db n for n = 1 . . . N. If these partial derivatives are all 
0 we have found the minimum, otherwise 

5. update a n and b n for n = l . . . N with 5a n and Sb n as in eqs. 26 and 27; 

6. restart from 2. 

This algorithm is known as steepest descent. In other approaches, like in the quasi-Newton 
methods, the local curvature of the functional is taken into account in order to compute an 
optimal descent step v and increase the convergence rate to the minimum. 

The adjoint-equations approach reduces the cost of computing the gradient; the cost is the 
same regardless of the number of controls. However, the number of flow-field evaluations to 
reach the minimum is very high, and grows more than linearly when the number of controls 
are increased, as shown by Beux & Dervieux [3]. Considering flows governed by the Navier- 
Stokes equations, as done by Cabuk h Modi [5] , it is seen that only a few minimization 
steps with a relatively coarse grid are possible, due to the excessive computational times 
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required. 

Ta’asan [25] proposed tlie pseudo-time method as a more efficient way of solving the op- 
timization problem. The main observation is the following. The solution of the optimization 
problem lies on the intersection of the hypersurfaces where the state, costate and optimality 
conditions are satisfied. Gradient based methods (including adjoint formulations) can be 
viewed as marching along the intersection of the state and costate hyper surf aces. This is 
an expensive process since at each step taken along the intersection requires the solution 
of two systems of PDEs. The idea of the pseudo-time method is to perform the marching 
on the optimality conditions hypersurface while relaxing the requirement to lie on the state 
and costate hypersurfaces. Note that the hypersurface representing the optimality condi- 
tions is one dimension less that of the state and costate equations. Therefore, the cost of 
such an iteration per step is significantly smaller than that of gradient based methods. Its 
convergence has been shown by Ta’asan [25] to be independent of the number of design 
variables. 

In some cases the optimality- conditions equations can be solved for the design variables 
and a simple implementation of the above idea exists. In other cases the optimality- 
conditions equations, viewed as an equation for the design variables keeping the state and 
costate fixed, may be singular and a more involved implementation is required. This is the 
case for the problem considered here. In such cases it is necessary to solve for the design 
variables together with the state and costate variables in a small neighborhood of the airfoil 
boundary, as shown by lollo, Kuruvila & Ta’asan [10]. The resulting algorithm is as follows: 
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1. start with a tentative set of a n and b n ] 

2. march the Euler equations, in time, a few steps on the entire field; 

3. march the adjoint equations, in time, a few steps on the entire field; 

4. solve in a small neighborhood of the airfoil boundary F the Euler equations and satisfy 
the boundary conditions, likewise do the same with the adjoint equations [10]; then 
compute the gradient of the functional; 

5. if the gradient is null, restart from step 2, repeating steps 3 to 5 until the state and 
costate equations are converged on the entire field. Otherwise update a n and b n with 
8a n and 8b n as in eqs. 26 and 27, go to 4. 

When this algorithm is used in conjunction with the inverse Theodorsen transform, con- 
vergence is obtained while increasing the number of coefficients used. Usually we started 
with N — 2 and increased N to the point where there is no appreciable decrease in the 
functional, usually for N = 6. All of the results presented in the next section are obtained 
by means of the pseudo-time method. The gradient is considered zero when it has decreased 
two orders of magnitude, while the functional S decreases differently according to the case 
considered. The computational cost of each optimization was about 5 times the cost of a 
single analysis of the flow field. 
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Figure 9: Pressure versus chord fraction for RAE2822 airfoil and optimal airfoil at M <*> = 
0.75. The unperturbed flow static pressure is taken as reference. 

6 Results and Discussion 

The analysis code was validated against known test cases. In particular we considered 
those contained in AGARD [1]. For example, we show the results for test case 02: a NACA 
0012 at 1 degree of incidence and = 0.85. The grid used had 240 tangential and 60 
radial points. The lift coefficient reported for this test case varies from 0.330 to 0.3889, the 
one that we find is 0.334. The drag coefficients reported range from 0.0464 to 0.0590; we 
have 0.0568. The iso-Mach lines of fig. 4 compare well with the benchmark of fig. 6.19 in 
AGARD [1]. 

In order to design an airfoil of some interest, it is necessary that it has some desirable 
aerodynamic characteristics, such as minimum drag, high lift, or given pitching moment. 
In addition some structural constraints, for example on the volume or on the trailing edge 
angle, must be satisfied. In principle, it would be possible to minimize wave drag for given 
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Figure 10: Pressure versus chord fraction for RAE2822 airfoil and optimal airfoil at M 0 0 = 
0.80. The unperturbed flow static pressure is taken as reference. 

lift using the optimization strategy that we have developed. The functional to be minimized 
would be a blend of drag and a quadratic penalty function for lift. 

We have preferred another approach that has the advantage of controlling the target 
pressure distribution produced on the airfoil. In fact, pressure distribution characteristics 
are important for boundary layer development and transition predictions, therefore the 
designer needs to control directly the pressure distribution over the airfoil. Furthermore, 
once pressure distribution is fixed, lift, drag, and pitching moment are loosely dependent 
on the actual geometry of the airfoil, such that if we define some parameters that control 
the distribution of pressure on the airfoil and we optimize such distribution for the desired 
aerodynamic characteristics, we have an a priori knowledge of the pressure target. 

Van Dam, van Egmond h Slooff [26], for example, propose that the velocity distribution 
over the upper and lower surface be divided in three regions: stagnation, followed by a rapid 
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Figure 11: Pressure versus chord fraction for RAE2822 airfoil and optimal airfoil at Moo = 
0.85. The unperturbed flow static pressure is taken as reference. 

acceleration; a region of slightly varying velocity that may be followed by a shock wave; and 
a pressure recovery region. These regions are defined by a number of fixed points in the 
pressure distribution, and the interpolation rules between these points take into account 
boundary layer development and separation requirements. 

In this work we have limited ourself to selecting as target pressure distributions of well 
studied airfoils, such as a Korn airfoil, the RAE 2822, and the NACA 0012. In particular, 
we have studied the problem of designing a family of wing sections that produce the same 
surface pressure, when varies. The same problem was studied for the first time by 
Chin [6] by means of inverse design, using the full potential equations. The airfoils he 
obtained had the surface pressure of a Korn profile at M <*> = 0.75. These airfoils had a 
trailing edge angle too small for practical applications, and in some cases the upper and 
the lower surfaces crossed each other. The next examples show the possibility of achieving 
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Figure 12: Supersonic bubble for AF-1 airfoil. 

given aerodynamic requirements, with practical airfoils. 

Using the same trailing edge of a NACA 0012 airfoil, we designed two airfoils at = 
0.75 and — 0.8 respectively, which approximate in a least square sense the pressure 
distribution of a Korn airfoil at — 0.7. See figs. 5-6. In fig. 7 are presented the airfoils 
obtained. In fig. 8 it is shown for comparison the pressure distribution over a Korn airfoil 
at Moo — 0.8. 

The maximum Mach number on the airfoil Mmax, the value of the functional at the 
minimum £ m j n5 the lift coefficient Cx, and wave drag Cd for each of the unperturbed flow 
conditions are reported in the next table. The drag coefficient was considered zero when 
it was of the same order of that corresponding to numerical dissipation for subsonic test 
cases, i.e. O(10"* 4 ). 
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Figure 13: RAE2822 airfoil and optimal designs at M ^ — 0.75, — 0.8 and M ^ = 0.85. 


AC, 

Afmax 

E • 
6 mm 

C L 

C D 

0.7 

0.98 


0.44 

0 

0.75 

0.99 

3.32 • 10- 4 

0.37 

0 

0.8 

1.04 

8.23 • 1CT 4 

0.31 

0 


Note that the difference in lift is a measure of the first order moment of the difference 
p(F) -p*. 

As another example of transonic design, the pressure generated over a RAE 2822 airfoil 
at Mqo = 0.7 was taken to be the target distribution. The trailing edge angle was required 
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to be the same as that of a NACA 0012, and three flight conditions were investigated: 
Moo = 0.75, Moo = 0.8 and = 0.85. The pressure distributions produced on the 
airfoils obtained approximated to different degrees the target pressure distribution. As 
Moo increases, the functional S has higher values at the minimum, i.e. for lower Moo the 
target and the obtained pressure distribution are closer, in a least square sense. Results 
for each case, see also figs. 9, 10, 11, 13, are shown in the next table. 


Moo 

Afmax 

< -'mm 

Cl 

Cd 

0.7 

0.98 


0.18 

0 

0.75 

1.00 

4.03 ■ 10~ 4 

0.30 

0 

0.8 

1.02 

5.86 ■ 10" 4 

0.26 

0 

0.85 

1.07 

1.03 • 10" 3 

0.22 

0 


The airfoil obtained for Moo = 0.85 is shown in fig. 11, and is designated AF-1. The large 
supersonic bubble generated on the profile is reported in fig. 12. The AF-1 airfoil has a 
maximum- thickness- to-chord ratio of 6.1% and, if compared to other transonic shock-free 
airfoils for = 0.85, it has favorable lift characteristics. In the work of Boerstoel [4], in 
fact, it is conjectured that the limit maximum-thickness-to-chord ratio for a shock-free non 
lifting airfoil it is about 7%. The nature of the recompression over the above airfoil was 
further investigated. The flow around the AF-1 airfoil was studied using a finer mesh with 
240 tangential and 60 radial points. The entropy rise across the weak shocks on upper and 
lower surfaces is 0(1 0 -3 ). The lift coefficient and the drag coefficient remain unchanged 
compared to the coarser-grid case. 
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It might be argued that the pressure distributions which were chosen for transonic design 
are rather peculiar: they show a typical plateau after an abrupt leading edge expansion 
combined with aft camber to get shockless recompression and for improving lift. For 
this reason, pressure distributions with other characteristics were tested, such as that of a 
NACA 0012 at 1 degree of incidence and Moo = 0.65. With this distribution, we designed 
at Mqo = 0.8 and the results are presented in figs. 14 and 15. We have: 


M x 

Mmax 

£ . 
nnin 

C L 

Cd 

0.65 

0.98 


0.29 

0 

0.8 

1.07 

3.79 • IQ" 3 

0.16 

0 


Compared to the other designs at M^ = 0.8 the approximation of the given pressure 
distribution is rather poor. 

From the results presented it is evident that the pressure distribution that can be approx- 
imated to a higher degree is that corresponding to the RAE2822 airfoil. 

7 Conclusions 

We have designed transonic wing sections of required performance using the Euler equa- 
tions. The wing sections designed have a given trailing edge angle and approximate, in a 
least square sense, a target pressure distribution. From a practical view point, the examples 
shown suggest that airfoils may be tailored for every transonic flight condition, taking into 
account constraints on the geometry. Furthermore we have shown that it is possible to 
generate profiles that approximate arbitrary pressure distributions, in particular it is pos- 
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Figure 14: Pressure versus chord fraction for NACA0012 airfoil at Moo = 0.65 1 degree of 
incidence and optimal airfoil at M^ = 0.65. 

sible to design profiles approximating a pressure distribution that satisfies requirements on 
pitching moment and boundary layer development. 

The design was conducted by means of an optimization algorithm that uses reasonable 
computational resources: one full optimization costs about 5 times the cost of one analysis. 
The algorithm is based on the pseudo-time method and on the inverse Theodorsen trans- 
form, that allows the minimization with an increasing number of harmonics, in a natural 
and consistent way. With this approach the design of aerodynamic components using the 
Navier-Stokes equations seems to be attainable. 
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Figure 15: NACA0012 airfoil and optimal design at = 0.8. 
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